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We describe the generalization of Analytic Perturbation Theory 
(APT) for QCD observables, initiated by Radyushkin, Krasnikov, 
Pivovarov, Shirkov and Solovtsov, to fractional powers of coupling 
— Fractional APT (FAPT). The basic aspects of FAPT is shortly 
summarized. We describe how to treat heavy-quark thresholds in 
FAPT and then show how to resum perturbative series in both 
the one-loop APT and FAPT. As an application we consider 
FAPT description of the Higgs boson decay — » bb. The main 
conclusion is: To achieve an accuracy of the order of 1% it is 
enough to take into account up to the third correction. 



r^, 1. APT and FAPT in QCD 



In the standard QCD Perturbation Theory (PT) we 
know that the Renormahzation Group (RG) equation 
das[L\/dL = —a^ — ... for the effective couphng 
a«(Q2) = a,[L]/Pf with L = HQ^A^), Pf = 
6o(iV/)/(47r) = (11 - 2iV//3)/(47rjl]. Then the one-loop 
solution generates Landau pole singularity, as[L] = 1/L. 

Strictly speaking the QCD Analytic Perturbation 
Theory (APT) was initiated by N. N. Bogolyuov et al. 
paper of 1959 [1], where ghost-free effective coupling for 
QED has been constructed. Then in 1982 Radyushkin 
[2] and Krasnikov and Pivovarov [3] using the same 
dispersion technique suggested regular (for s > A^) QCD 
running coupling in Minkowskian region, the famous 

1 TT 

— arctan— . After that in 1995 Jones and Solovtsov 

TT L 

discovered the coupling which appears to be finite for 
all s and coincides with Radyushkin one for s > A^, 
namely 2li [L] in Eq. (|2b|) . Just in the same time Beneke 
et al. [4,5], using the renormalon based approach, and 
Shirkov and Solovtsov [6], using the same dispersion 
approach of [1], discovered ghost-free couphng Ai [L] , Eq. 
l(2a|) . in Euchdean region. 

But Shirkov-Solovtsov approach, named APT, was 
more powerful: in Euclidean domain, —q^ = Q^, L = 
InQ^/A^, it generates the following set of images for 
the effective coupling and its n-th powers, {^„[L]}^gj^, 
whereas in Minkowskian domain, = s, Lg — Ins/A^, 
it generates another set, {2tn[is]}„gp^ (see also in [7]). 



APT is based on the RG and causality that guaranties 
standard perturbative UV asymptotics and spectral 
properties. Power series X^m'^'narl^] transforms into 
non-power series dmAm[L\ in APT. 

By the analytization in APT for an observable f{Q^) 
we mean the "Kallen-Lehman" representation 
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with Pf{<y) = — Im [/(— <t)] . Then in the one-loop 



approximation pi (cr) 
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whereas analytic images of the higher powers {n> 2, n e 
N) are: 
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d_y-^ Ai[L\ 
(n- 1)! V dL) \^i[Ls\ 

glance, the APT is a complete theory 
to produce an analytic answer for 
perturbative series in QCD. But in 2001 
Karanikas and Stefanis [8] suggested the principle of 
analytization "as a whole" in the plane for hadronic 
observables, calculated perturbatively. More precisely, 
they proposed the analytization recipe for terms hke 
J^dxJ^dy as (Q^xy) f{x)f{y), which can be treated as 
an effective account for the logarithmic terms in the 
next-to-leading-order approximation of the perturbative 
QCD. This actually generalizes the analytic approach 
suggested in [9]. Indeed, in the standard QCD PT one 
has also: 

(i) the factorization procedure in QCD that gives rise 
to the appearance of logarithmic factors of the type: 
<[L]L; 

(ii) the RG evolution that generates evolution factors of 
the type: B{Q^) = [Z{Q^)/Z{ii^)] B{fi^), which reduce 



in the one-loop approximation to Z{Q^) ^ a^lL] with 

f{Q^) and f[L] in order to specify the arguments we mean — squared momentum or its logarithm 
ln(Q2/A2), that is f[L] = /(A^ ■ e^) and A^ is usually referred to Nf = 3 region. 



We use notations 
— l„i'r)2 /a2 
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V = 7o/(26o) being a fractional number. 
All that means that in order to generalize APT 
in the "analytization as a whole" direction one 
needs to construct analytic images of new functions: 
a^, agL"*,.... This task has been performed in the 
frames of the so-called FAPT, suggested in [10,11]. Now 
we briefly describe this approach. 

In the one-loop approximation using recursive 
relation ([3]) we can obtain explicit expressions for Av \V\ 
and 21^ [L]: 



F{e 
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Here F(z, v) is reduced Lerch transcendental function, 
which is an analytic function in v. They have very 
interesting properties, which we discussed extensively in 
our previous papers [10-13]. 

Construction of FAPT with fixed number of quark 
flavors, iV/, is a two-step procedure: we start with 
the perturbative result [as((5^)] , generate the spectral 
density Pi/(ct) using Eq. ([l]), and then obtain analytic 
couplings Au[L] and 2ty[i] via Eqs. Here Nj is 
fixed and factorized out. We can proceed in the same 
manner for 7V/-dependent quantities: [a^lQ'^; Nf)'^ ^ 
p^(cr;iV/) = pALa^Nf] = pA<y)/f3f A[L;Nf] and 
2lj^[L; Nf] — here Nf is fixed, but not factorized out. 

Global version of FAPT [12], which takes into 
account heavy-quark thresholds, is constructed along 
the same lines but starting from global perturbative 
coupling [a|'°''((52)j ^ being a continuous function of 
Q2 due to choosing different values of QCD scales Af, 
corresponding to different values of Nf. We illustrate 
here the case of only one heavy-quark threshold at s = 
ml, corresponding to the transition Nf = 3 ^ Nf — i. 
Then we obtain the discontinuous spectral density 



(a) =9{L^<L4) pn[Ly,3] 

+ e{Li< La) Pn [La+Xi-A] , 



(5) 

with La = ln(o-/A§), Lf = Inlmj/Alj and A/ = 

ln(A2/A2) for / = 4, which is expressed in terms 

of fixed-fiavor spectral densities with 3 and 4 fiavors, 
pn[L;3] and pn[L + A4;4]. However it generates the 
continuous Minkowskian coupling 

Slf'lL] ^e{L<L4) (2l,[L;3]-h A43a,) 
+ 9{L4<L) 2l^[L-f A4;4]. 



with A432l^ = 2l^,[X4 + A4; 4] -2l^,[L4; 3] and the analytic 
Euclidean coupling Af!°'' [L] 



Af;°AL] =^[L-f A4;4] 

Li _ 

p^La] 3] - p^La + A4;4] 



dL„ 



(6b) 
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(for more detail see in [12]). 

2. Resummation in the one-loop APT and 
FAPT 

We consider now the perturbative expansion of a typical 
physical quantity, like the Adler function and the ratio 
R, in the one-loop APT. Due to limited space of our 
presentation we provide all formulas only for quantities 
in Minkowski region: 



n[L] = Y,dn%^[L] 



(7) 



We suggest that there exist the generating function P{t) 
for coefficients dn = dn/di: 



d„ = / P{t)e-^dt with 
'o 



P{t) dt = l. 



(8) 



To shorten our formulae, we use for the integral 
J^f{t)P{t)dt the following notation: ((/(t)))p(t). Then 
coefficients c?„ — di ((t"^^))p(t) and as has been shown 
in [14] we have the exact result for the sum in ([7]) 



7^[L] = dl((2ll[L-^]))p(,) 



(9) 



The integral in variable t here has a rigorous meaning, 
ensured by the finiteness of the coupling 2li [t] < 1 and 
fast fall-off of the generating function P{t). 

In our previous publications [12, 15] we have 
constructed generaHzations of these results, first, to the 
case of the global APT, when heavy-quark thresholds 
are taken into account. Then one starts with the 
series of the type lO, where 2t„[L] are substituted by 
their global analogs 21^'°'' [L] (note that due to different 
normalizations of global couplings, 21^'°'' [L] ~ 2l„[L]//?/, 
the coefficients c?„ should be also changed). Then 



TZ'""'[L] = di 6i(i < L4) ( (A42I1 M -fSti 



))p(t) 



-die{L>Li){{^] 



L + M-—-A 

PA 



))p[t) ; (10) 



(6a) where A42l^[i] =21^ L4 + A4-i//34;4 -21^ La-i/Z^a; 3 . 
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The second generalization has been obtained for the 
case of the global FAPT. Then the starting point is 
the series of the type X^^^o ^" ^n+i'I^] result 
of summation is a complete analog of Eq. (flOl) with 
substitutions 



1 



(11) 



do ^ do 21,. [L] ,^i[L-t]^^i+^[L-t], and A42I1 [t] ^ 
A42li+^[t]. All needed formulas have been also obtained 
in parallel for the Euclidean case. 

3. Applications to Higgs boson decay 

Here we analyze the Higgs boson decay to a bb pair. For 
its width we have 



r(H ^ 66) 



MhRs{MI) 



(12) 



with Rs{Mjj) = mliMfj) Rs{Mfj) and i?s(s) is the R- 
ratio for the scalar correlator, see for details in [10, 16]. 
In the one-loop FAPT this generates the following non- 
power expansiorll: 



and 2lf°''(TO^) = 0.122. In this range we have (Lg = 
In(m2/Ai)) 



3 rh 



(1) 
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(15) 



with P^^it) defined via Eqs. dH]) and (fTTj) . 

Now we analyze the accuracy of the truncated FAPT 
expressions 



TZs[L;N] = 3i 



^ Js 



(16) 



and compare them with the total sum TZs [L] in Eq. (fTSl 
using relative errors AAr[L] = 1 — TZs[L; N]/TZs[L]. In 
Fig. 1 we show these errors for = 2, iV = 3, and 
= 4 in the analyzed range of L G [11, 13.8]. We see 
that already TZg [L; 2] gives accuracy of the order of 2.5%, 
whereas TZs[L; 3] of the order of 1%. 



7^s[L] = 3i 



*(1) 



[L] 



E 



(13) 



where m^^^ 



9.05 ± 0.09 GeV^ is the RG-invariant of 
the one-loop rnKfj,"^) evolution ml{Q^) = m^-^^ a^" (Q^) 
with = 270/60(5) = 1.04 and 70 is the quark-mass 
anomalous dimension. This value m^^^ has been obtained 
using the one-loop relation [17] between the pole 6-quark 
mass of [18] and the mass mb(mb). 

We take for the generating function P{t) the Lipatov- 
like model of [15] with {c = 2.4, f3 = -0.52} 



Ps{t) 



jt/c) + (3 



i) + /?r(n) 



1 + p 

-t/c 



(14a) 
(14b) 



It gives a very good prediction for with n = 2,3,4, 
calculated in the QCD PT [16]: 7.50, 61.1, and 625 in 
comparison with 7.42, 62.3, and 620. Then we apply 
FAPT resummation technique to estimate how good 
is FAPT in approximating the whole sum 7?.s[L] in 
the range L g [11.5, 13.7] which corresponds to the 
range Mh £ [60, 180] GeV^ with Aq^^^ = 189 MeV 




12 12.5 13 13.5 

Fig. 1. The relative errors Ajv[L], A'' = 2,3 and 4, of the truncated 
FAPT in comparison with the exact summation result, Eq. ||15|I . 

Looking in Fig. 1 we understand that only in order 
to have the accuracy better than 0.5% we need to take 
into account the 4-th correction. We verified also that 
the uncertainty due to P(<)-modelling is small < 0.6%, 
while the on-shell mass uncertainty is of the order of 2%. 
The overall uncertainty then is of the order of 3%, see 
in Fig. 2. 



^Appearance of denominators tt" in association with the coefficients dn is due to dn normalization. 
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[ Mh [GeV] : 

100 120 140 160 180 

Fig. 2. The width Fji^ as a function of the Higgs boson mass 
Mh in the resummed FAPT. The width of the shaded strip is due 
to the overall uncertainties, induced by the uncertainties of the 
resummation procedure and the pole mass error-bars. 

4. Conclusions 

In this report we described the resummation approach 
in the global versions of the one-loop APT and FAPT 
and argued that it produces finite answers, provided the 
generating function P{t) of perturbative coefficients d„ 
is known. The main conclusion is: To achieve an accuracy 
of the order of 1% it is enough to take into account 
up to the third correction— in complete agreement with 
Kataev&Kim [17]. The di coefficient value is needed only 
to estimate corresponding generating functions P{t). 
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